Dissipation induced Instabilities and the Mechanical Laser 
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We study the 1:1 resonance for perturbed Hamiltonian systems with small dissipative and energy 
injection terms. These perturbations of the 1:1 resonance exhibit dissipation induced instabilities. 
This mechanism allow us to show that a slightly pumping optical cavity is unstable when one takes 
into account the dissipative effects. The Maxwell-Bloch equations are the asymptotic normal form 
that describe this instability when energy is injected through forcing at zero frequency. We display 
a simple mechanical system, close to the 1:1 resonance, which is a mechanical analog of the Laser. 
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A fundamental problem in mechanics is the determi- 
nation of the stability of equilibria. Hamiltonian vector 
fields can undergo a variety of instabilities as a single 
bifurcation parameter is varied Q] . There are two funda- 
mental codimension-one bifurcations of Hamiltonian sys- 
tems. The first is the stationary or steady state bifur- 
cation, which is characterized by two eigenvalues merg- 
ing at zero with multiplicity two Q. The second is the 
1:1 resonance, which is the collision of two pure imagi- 
nary eigenvalues (and their complex conjugates) at finite 
frequencies with multiplicity two Q. This latter insta- 
bility has been called by different names depending on 
the field where it has appeared; for instance, 'confusion 
of frequencies' in mechanical or electrical oscillators || , 
'dispersive instability' Q and 'Hamiltonian Krein-Hopf 
bifurcation' ||. These instabilities are a consequence of 
the fact that in the Hamiltonian case, if A is an eigen- 
value, then H so is —A. The same property is present in 
time reversible systems, i.e., they exhibit identical generic 
instabilities. Recently, the instabilities of quasi- reversible 
systems have also been characterized J^, |t) , in which the 
irreversible effects are small and can be considered as 
perturbative terms close to the instability. In particular, 
systems close to the quasi-reversible 1 : 1 resonance are de- 
scribed by the Maxwell-Bloch equations when the energy 
is injected through a forcing at zero frequency (?]]. The 
Maxwell-Bloch equations describe the interaction of an 
electromagnetic field and a collection of two level atoms 
at an optical cavity 

The aim of this letter is to study the 1:1 resonance for 
systems in the neighborhood of a Hamiltonian one; that 
is, we shall consider perturbed Hamiltonian systems with 
dissipative and energy injection terms. Near this insta- 
bility, the dissipative terms are responsible for a spectral 
bifurcation, i.e., the dissipation induces an instability 
This mechanism allows us to show that a slightly pump- 
ing optical cavity is unstable when one takes into account 
the dissipative effects. The Maxwell-Bloch equations are 
the asymptotic normal form that describes this insta- 
bility in the presence of a conserved quantity. We shall 
display a simple mechanical system, which we call the 
Mechanical laser, which, close to the 1:1 resonance, is a 
mechanical analog of the Laser. 

For the sake of simplicity, we first consider a Hamilto- 



nian system with two degrees of freedom H [p l , qi, {A}), 
where p l and qi are canonically conjugate variables, 
(i = 1,2), and {A} is a set of parameters. Assume there 
is a 1:1 resonance for an equilibrium at A = A c , i.e., the 
spectrum at the equilibrium has a pair of pure imaginary 
eigenvalues of multiplicity two, say at ±iQ. Nearby, the 
instability the system is governed by the normal form [Tc| ] 

d t A = iQA + B 

d t B = mB + f(\A\ 2 ,i(AB*-BA*),{X-X c })A + 
ig(\A\ 2 ,i{AB* -BA*),{\-\ C })B, (1) 

where / and g are complex polynomial functions. There 
is a change of variables from the given ones {p 1 ,qi,p 2 ,q2j 
to the new ones (the complex variables A, B) of the form 
{p 1 , qi , p 2 , q 2 } = A$ + Bx + N (A, B, A* ,B*), where * is 
the eigenvector of the linearized system corresponding to 
iQ and x is a generalized eigenvector in the Jordan sense, 
and N (A, B, A*, B*) is a nonlinear vector. This change 
of variables is not canonical in general. When one con- 
siders rotated variables [A — e int A',B = e l0t <9 t A') and 
the dominant terms, the normal form reads (omitting 
the primes) 



d tt A = eA + idd t A - a\A\ 2 A 



(2) 



where e is the bifurcation parameter, which is propor- 
tional to A — A c ; henceforth we assume e <C 1, S is the 
gyroscopic term jj) or detuning term || and a is an or- 
der one parameter. The variables and parameters scale 
as A ~ e 1 ' 2 , dtA ~ e, d tt A ~ e 3 / 2 , and 8 ~ e 1 ' 2 . All 
terms in this equation are of order e 3 / 2 , while the higher 
order terms are e 2 . 

When equation (|l|) is written as a second order system, 
it has an imaginary term linear in A, but such a term is 
not present in (j^) because it breaks the eigenvalue sym- 
metry (A — > —A). The asymptotic normal form (^) has 
the following Hamiltonian 



H = d t Ad t A* - e\A\< 



(3) 



with the Poisson-Bracket (F and G are real) 
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Thus, to cubic order, the nonlinear change of variables is 
canonical. The eigenvalues of the zero solution (A = 0) 
are ±i (8/2 ± 1/2 V<5 2 - 4e). When £ - (5 2 /4 is nega- 
tive, the initial Hamiltonian system has four distinct pure 
imaginary eigenvalues; e — J 2 /4 equal to zero is the 1:1 
resonance at ±f2 frequencies and when it is positive, the 
eigenvalues have nonzero real part. Note that the gyro- 
scopic term is a stabilizing effect ||. 

We now consider this Hamiltonian system under the 
influence of small dissipative terms. This leads to a new 
term in the asymptotic normal form as follows: 

d t tA = eA-(fi-i6)d t A-a\A\ 2 A, (4) 

where fx is positive and order e 1 / 2 . 

To study the effects of dissipative terms in ([| ) , we con- 
sider its characteristic polynomial, which has the form 
A 4 + 2/xA 3 + (2e - y 2 + 5 2 ) X 2 - A2e/i + e 2 . Looking for 
roots of the form A = a±ib and A = c±id, one recognizes 
that (J, — — (a + c) and Efi = (c (a 2 + b 2 ) + a (c 2 + d 2 )). 
Hence, when e is negative, a and c are negative, i.e., all 
eigenvalues are to the left of the imaginary axis. For e 
positive and e — 5 2 /4 negative, the unperturbed system 
is marginal, but the perturbed one satisfies ac < 0; in ad- 
dition, the eigenvalues with larger frequency move to the 
left of the imaginary axis (stable modes) and the others 
to the right (unstable modes), but the eigenvalues that 
move furthest away from the axis are stable. Finally, 
when e > S 2 /A, the eigenvalues have nonzero real part, 
and again the stable modes are the furthest from the 
imaginary axis. Therefore, we are led to consider dissi- 
pation induced instabilities ||. The destabilizing effects 
through positive or negative total dissipative perturba- 
tion was known a long time ago by Lord Kelvin et al. || 
From the physical point of view, one can understand this 
phenomena as follows: when e is negative, the energy (||) 
has a minimum at the origin, hence when dissipation is 
added, the solutions near the origin move towards it. By 
contrast, if e is positive, the energy has a saddle point at 
the origin and when e < S 2 /A, this is unstable with an 
algebraic evolution in time, so when one adds dissipation, 
this is consistent with the conclusion that the solutions 
near the origin move away from it, exponentially in time. 

Using the preceding analysis, we infer that close to the 
1:1 resonance, generically the dissipative terms induce 
an instability. In the case that the instability happens 
with null detuning (6 = 0), the normal form has only a 
real coefficient and so the dissipative terms do not induce 
instability. A physical example of this last situation is the 
weakly dissipative baroclinic instability when the effect 
of earth's sphericity is ignored Jll| and another is the 
Kelvin-Helmholz instability |fL2| . 

As we shall see in detail later, the Laser is a system 
that shows dissipation induced instability, but firstly, we 
need to discuss how the energy is injected to the system, 
so that the solution that becomes unstable is persistent 



when the dissipative terms are added. There are two nat- 
ural ways to inject energy to the modes, namely through 
forcing at finite frequency or at zero. The latter situation 
is common in physical systems. In a Hamiltonian system 
this is only possible if there is a conserved quantity, that 
is, a zero eigenvalue whose mode is nonlinearly coupled 
with the other ones. For instance, when there is a cyclic 
variable, the respective momentum Z is conserved. A 
Hamiltonian system that has a 1:1 resonance in the pres- 
ence of conserved quantity, leads us to consider an extra 
equation of the form dtZ — and an additional term 
— ZA in the equation ( Q). The system presents different 
behaviors depending on the value of e — Z . When one 
includes dissipative terms and forcing, the asymptotic 
normal form reads 

d tt A = (e- Z)A-{fi-i8)d t A~a\A\ 2 A 
d t Z = vZ + i-,\A\ 2 (5) 

where Z ~ e, v ~ e 1 ! 2 and r\ ~ e 1 / 2 . The term r]\A\ 2 per- 
mits a non trivial coupling between the variables. When 
the unperturbed Hamiltonian system has more modes 
without resonances between them, the perturbed system 
is governed by the above equations, since the intensi- 
ties of the other modes decreases in time. Through a 
nonlinear change of variables, the previous equations are 
equivalent to the Maxwell-Bloch equations 0. Using a 
multiscaling method, the dispersive instability with small 
dissipation is also described by the previous equations f| . 

To illustrate how dissipation induced instabilities en- 
ter, we consider the semiclassical description of the Laser. 
This is based on the self-consistent interaction of the 
electromagnetic field with an active medium within an 
optical cavity. The electric field is described classically 
(by the Maxwell equations) and the matter as ensem- 
ble of atoms possessing two quantized energy levels; phe- 
nomenological terms are added to complete the descrip- 
tion. Thus, the system is described by |j, [ll| 

d 2 E _ d 2 E d 2 P dE 

^- = -^-hi + (i + s) 2 )?-9NE, 

dN fdP \ 

— =- lll (N- No )+E^—+ 1± Py (6) 

with periodic boundary condition at the cavity length (L, 
E(0,t) = E(L,t)). Here E, P and N are dimensionless 
quantities, which correspond to linearly polarized elec- 
tric field, the dipole polarization field and the population 
inversion. 711 , 7j_ are the decay rate associated to sponta- 
neous emission and interaction between the atoms, K is a 
damping related to the mirror losses, S is the detuning, g 
is a coupling constant which characterizes the atoms and 
Nq the pump parameter. In the time reversible limit of 
the above equations, i.e., 7^ = 7|| = n = 0, the system 
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has Hamiltonian density |b 

fD 1 , 
n=(- Re (a+a*_) 

with the Poisson-bracket 
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Where N = \a + \ 2 - |a_| 2 , P = ^Re(a+a*), £ = A t , 
and w = I + 6. The Hamiltonian is just the sum of the 
electromagnetic energy and the atomic excitation energy 

B 

Changing the cavity length in the time reversible limit 
of equation (^), leads to a 1:1 resonance for the non- 
lasing solution (E = P = 0, N = Do), which gives rise 
to an electromagnetic wave with ±oj frequencies. Using 
the slowly varying envelope (WKB) approximation leads 
to the Maxwell-Bloch equations 0] . 

Figure 1 shows the space-time diagram of the electric 
field of a numerical simulation of the semiclassical model, 
plus noise with small intensity, close to the 1:1 resonance 
without dissipative terms (see fig. la) and with dissi- 
pative ones (see Fig. lb). The numerical simulations 
start with the same initial condition, namely, the no- 
lasing solution with excited atoms (D > 0). It is clear 
from these pictures that the inclusion of dissipation in- 
duces the Laser to respond. One can physically under- 
stand what happens, since without dissipation the atoms 
excited decay through the stimulated emission, i.e., the 
non-lasing solution becomes unstable very slowly (non- 
linear mechanism). Instead, when one takes into account 
the dissipative terms, the excited atoms decay, for in- 
stance through stimulated emission and collisions, i.e., 
exponentially in time. Note that, one must to pump to 
No = D in order that the non-lasing solution persists. 
In brief, with nonzero detuning and slight pumping, the 
dissipation induces the Laser to respond. 

To illustrate the 1:1 resonance in a simple Hamiltonian 
system, we consider a mechanical system, that we call the 
Mechanical Laser, which consists of two coupled spherical 
pcndula in a gravitational field, with a support, which 
can rotate around a vertical axis. The lower pendulum is 
constrained to move in a plane that is orthogonal to the 
plane of the upper pendulum (see Fig. 2). 

The system rotates with angular velocity <p with re- 
spect to the vertical. The quantities mi, m 2 , h and 1% 
are the mass and length of the upper and lower pendula 
respectively, and I is the dimensionless moment of inertia 
of the support. The system will dissipate energy because 
of friction at the contacts and the motion of the pendu- 
lum masses in a fluid (for example, the air) via Stokes' 
law. Energy is injected through a constant torque at the 
upper pendulum pivot point. The governing equations 



for the angles Q\ ft) and #2 (t) and the vertical angular 
velocity tp read [|l5) 

#1 = — a 2 sin 9\ sin — c 2 sin 9\ cos #2# 2 

—2a 2 cos 9\ cos6*2^#2 + sin^i cos^i^j 2 

2 9 n 

—a" cos 9\ sin 02<p — — sin 6\ — v\9\, 

62 = — sin 6*i sin 6261 — cos 9\ sin 62O 2 
+2 cos Q\ cos 02(p&i + sin #2 cos 02'p 2 
+ sin0i cos # 2 <^ — y sin 6*2 — ^2^2, 



d 
dt 



(sin 



2 " 1 a 2 sm 2 2 )<p- 



—a 2 sin 61 cos 62Q2 + I<p 



a 2 cos^i sin 82O1 



= —v^ {tp — ft) — fix sin 2 d\<p — \i2 (sin 2 9\ + sin 2 62) tp. 

where V\, V2, v v , Ml an d M2 are damping coefficients, 
a = \J TO2/ (m~i + 7712) is the relative factor of the energy 
between the oscillators, and we have written the torque 
as u v fl. For the sake of simplicity, we have considered the 
case of pendula of equal lengths (li = fe = I). When one 
regards the Hamiltonian limit of the previous equations, 
the vertical solution or non-lasing solution 9i =62 = 0, 
ip = Qq has a 1:1 resonance when 



^0 — — 1/ T 




with frequency uj c = ±y/gm2/lmi. The centripetal force 
is more intense than the gravitational force when f^o > 
f2 5 = y/g/l. As a consequence, the Coriolis force exerted 
by one pendulum on the other, the non-lasing solution 
is marginal when Q g < £lo < Q c . In this region the 
system is nonlinearly unstable and the system becomes 
linearly unstable when £Iq > £l c ; this exhibits a coherence 
oscillation, which is the signatures of the laser instability. 

Near the 1:1 resonance, the coupled pendulum is gov- 
ern by (H), where fill 



o(0-Q c ) g /a 4 - 2cr 3 - 2cr 2 + 3 

l VL r ' 42 V 1 - a 2 



S 
V 



2a{n-fl c ), v 



Mi 



1 

21 

v in ( SX 



2 {"I + V2) ■ 



1 + 32)M 2 -2^I^-1 



where the variables are related to the dominate order by 
A = (ct0i + id 2 )evp(iuj c t)/2 and Z = 2VL c o 2 (fi c - tp). 
Transforming equations fJJ) to the Maxwell-Bloch equa- 
tions, one obtains the following analogy of the electric 



field, polarization and population inversion 



E = 

P = 

N = 



{ct#i + ie 2 } 



e t(A+u e )t + 



{<j6 x + i6 2 +a(a9 1 +i6 2 )\ 



(A+uj c )t 



+ C.C. 



2fl r a 2 



a 



a 2 e\ + el 



where A = 5 (ya + rf) / (2va + 2rj — fia) , a ~ 
(va + rj) I a + i2A and e = (e + A 2 - 6 A) ja + 
[ya + r\) (fia — va — r])/a 3 . The mechanical analog of 
the electric field is a simple function of the pendulum 
displacements with respect to the vertical. In the Hamil- 
tonian case the polarization is the slow time derivative of 
the electric field. The mechanical analogy of the stimu- 
lated emission mechanism (time reversal effect) is the ver- 
tical angular momentum conservation. Thus, when one 
increases the pendula tilt (the intensity of the electric 
field enlargement), the angular velocity decreases (the 
population inversion declines). This is a main ingredi- 
ent of Laser theory. If the torque or pumping is bigger 
than the gravitational force (f2o > ^g)> i- e -j the system is 
excited in Laser terminology, then the torque that gives 
rise to an angular velocity slightly higher than fl g leads 
to a dissipation induced instability. 

In summary, we have studied the 1:1 resonance for per- 
turbed Hamiltonian systems with small dissipative and 
energy injection terms. Nearby, the 1:1 resonance ex- 
hibits dissipation induces instabilities. This allow us to 
show that a slightly pumping optical cavity is unsta- 
ble when one takes into account the dissipative effects. 
The Maxwell-Bloch equations are the asymptotic normal 
form that describe this instability when energy is injected 
through forcing at zero frequency. We have displayed a 
simple mechanical system, the Mechanical laser or dou- 
ble spherical pendulum with a support, which, close to 
the 1:1 resonance, is a mechanical analog of the Laser. 
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FIG. 1: Spatio-temporal diagram of the electric field of the semi- 
classical model (equation (Q)) with L = 120, 8 = 0.94, g = 0.6, 
D = 0.0338 a) without dissipation n = -fj_ = 711 = 0, b) with 
dissipation re = 0.03, 7 ± = 0.01, 7y = 0.03 and N a = 0.0338. 




FIG. 2: Schematic representation of the Mechanical Laser. 



